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, Abstract. The aim of this paper is to contribute to understanding of the pattern for- 

' mation phenomenon in reaction-diffusion equations coupled with ordinary differential 

s ! , equations. Such systems of equations arise, for example, from the modeling of intcr- 

1 actions between proliferating cells and diffusing growth factors. We focus on stability 

& , of solutions to an initial-boundary value problem for a system consisting of a single 

reaction-diffusion equation coupled with an ordinary differential equation. We show that 
0^ | such systems exhibit diffusion-driven instabilities (Turing instability) under the condition 

of autocatalysis of non-diffusing component. Nevertheless, there exist no stable Turing 
p I , patterns, i.e. all continuous spatially heterogenous stationary solutions are unstable. In 

■ addition, we formulate instability conditions for discontinuous patterns for a class of 

nonlinearities. 

> 

a ■ 
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1 — 1. Introduction 

Diffusion- driven instability (DDI), called also Turing-type instability, is a phenomenon 
CN| . in mathematical biology, which has been often used to explain de novo pattern forma- 
tion. Let us recall that DDI is a bifurcation that arises in a reaction-diffusion system, 
| when there exists a spatially homogeneous solution, which is asymptotically stable with 
respect to spatially homogeneous perturbations, but unstable to spatially heterogeneous 
perturbations. Models with DDI describe then a process of a destabilization of stationary 
spatially homogeneous steady states and evolution of spatially heterogeneous structures 
towards spatially heterogeneous steady states. The ideas on DDI have inspired develop- 
' ment of a vast number of mathematical models since the seminal paper of Turing [36], 
providing some explanations on symmetry breaking and de novo pattern formation dur- 
^ ! ing development, explaining shape of animal coat markings, and predicting oscillating 
chemical reactions, see e.g. [22] and references therein. 

In particular, the following system of reaction-diffusion equations 

u t = eAu + f(u, v), v t = DAv + g(u, v), 

considered in a bounded domain and supplemented with boundary and initial conditions, 
has been proposed as a mathematical model describing the diffusion-driven instability. 
Here, the unknown function u = u(x,t), sometimes called the activator, diffuses with a 
very small diffusion coefficient e > 0. On the the hand, v = v(x, t) describes the density of 
the inhibitor and the corresponding equation contains a large constant diffusion coefficient 
D > 0. 

In several models, however, the diffusion of u is either so small that it is practically 
negligible or u does not diffuse at all, which leads to a system of one ordinary differential 
equation coupled with one reaction-diffusion equation. Such systems of equations arise e.g. 
from modeling of interactions between cellular processes and diffusing growth factors, see 
[18, 13] for derivation of such mathematical models, [9, 12, 15, 16, 17] for some biological 
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applications and [14] for their mathematical analysis. Such models are different from 
classical Turing-type models and the spatial structure of the pattern emerging from the 
destabilization of the spatially homogeneous steady state cannot be determined based on 
linear stability analysis. Therefore, the existence and stability of spatially heterogeneous 
patterns arising in models exhibiting diffusion- driven instability, but consisting of only 
one reaction-diffusion equation is an interesting issue. As shown in [14] for a particular 
example of a react ion-diffusion-ODE system, it may happen that there exist no stable 
Turing-type patterns and the emerging spatially heterogenous structures have a dynamical 
character. In numerical simulations, solutions having the form of periodic or irregular 
spikes have been observed. 

The aim of this work is to investigate to which extent these results, obtained in [14], con- 
cerning the instability of all stationary structures can be generalized to react ion-diffusion- 
ODE models with a single diffusion operator. 

In the remainder of this paper, we focus on the following two-equation system 

(1.1) u t = f(u,v), for x e fi, t > 0, 

(1.2) v t = Av + g(u,v) for x e Q, t>0 

in a bounded domain Q C WL N for iV > 1, with a sufficiently regular boundary dQ, 
supplemented with the Neumann boundary condition 

(1.3) d v v = for xeffi, t> 0, 

where d u = djdv and v denotes the unit outer normal vector to dQ, and initial data 

(1.4) u(x , 0) = u (x) , v(x, 0) = Vq(x). 

The nonlinearities / = f(u,v) and g = g(u,v) are arbitrary C 2 -functions that satisfy 
certain natural (biologically motivated) assumptions. Notice also that the diffusion term 
in equation (1.2) could have the form DAv with a constant diffusion coefficient D > 0, 
however, after rescaling the time variable and the functions / and g, we can assume that 
D = 1. 

By a standard theory (see e.g. [29] and the beginning of Section 4), the boundary value 
problem (1.1)-(1.4) has a unique local-in-time solution e.g. for every uq, vq 6 L°°(Q). The 
aim of this work is to investigate the stability properties of stationary solutions. 

The paper is organized as follows. In Section 2, main results of this paper are formu- 
lated. Section 3 provides examples of reaction-diffusion-ODE systems considered in math- 
ematical biology, which fit the framework of the theory developed in this work. Proofs 
of the results are presented in Sections 4 and 5. Section 4 is devoted to showing insta- 
bility of continuous stationary solutions under so-called autocatalysis and compensation 
conditions. Moreover, for the instability of discontinuous solutions additional conditions 
on the model nonlinearities are required. In Section 5, constant and non-constant sta- 
tionary solutions of considered models are characterized and it is shown that they satisfy 
the conditions formulated in Section 4. Appendix contains additional information on the 
model of early carcinogenesis which was the main motivation for the research reported in 
this work. 

2. Results and comments 



2.1. Instability of steady states. The main goal of this paper is to formulate assump- 
tions on stationary solutions of problem (1.1)-(1.4), which lead to their instability. Then, 



UNSTABLE PATTERNS IN REACTION-DIFFUSION-ODE SYSTEMS 



3 



we show that such assumptions are satisfied in a natural way in reaction-diffusion-ODE 
systems modeling the Turing instability. 

First, we consider regular stationary solutions (U, V) of problem (1.1)-(1.3), namely, we 
assume that there exists a solution (not necessarily unique) of the equation f(U (x), V(x)) = 
that is given by the relation U{x) = k(V(x)) for all x G Q with a C 1 -function k = k(V). 
Thus, every regular stationary solution (U, V) of the boundary value problem 



(2.1) f(U,V) = for xgH, 

(2.2) AV + g(U,V) = for ifll, 

(2.3) d u V = for x G d£l 
satisfies the elliptic problem 

(2.4) AV + h(V) = for xefi, 

(2.5) d u V = for x G dSl, 
where 

(2.6) h(V) = g(k{V),V) and U(x) = k(V(x)). 



Regular stationary solutions appear to be unstable solutions to problem (1.1)- (1.4) 
under a simple assumption imposed on the first equation. 

Theorem 2.1 (Instability of regular solutions). Let (U, V) be a regular solution of problem 
(2.1)-(2.3) satisfying the following "autocatalysis condition": 

(2.7) f u {U(x), V(x)) > for all x G H. 

Then, (U,V) is an unstable solution the initial-boundary value problem (1.1)-(1.4). 

Inequality (2.7) can be interpreted as an autocatalysis in the dynamics of u at the 
steady state (U,V). Notice that, in this work, the stability of a stationary solution 
is understood in the Lyapunov sense. Moreover, we prove the nonlinear instability of 
stationary solutions to problem (1.1)-(1.3) and not only their linear instability, i.e. the 
instability of zero solution of the corresponding linearized problem, see Section 4 for more 
explanations. 

Each constant solution (u,v) G 1R 2 of problem (2.1)-(2.3) is a particular case of regular 
solutions. Thus, Theorem 2.1 provides a simple criterion for the diffusion-driven instability 
of (u, v). 

Corollary 2.2. If a constant solution (u,v) of problem (1.1)-(1.4) (namely, f(u,v) = 
and g(u,v) = 0) satisfies the inequalities 

(2.8) Uu, v) > 0, Uu, v) + gJu, v) < 0, det ( V \ ^ V \ ] > 0, 

then it has the DDI property. 

This corollary results immediately from Theorem 2.1, because the second and the third 
inequality in (2.8) imply that (u,v) is stable under homogeneous perturbations, see Re- 
mark 5.2 below for more details. 

Remark 2.3. The instability results from Theorem 2.1 and Corollary 2.2 can be summa- 
rized in the following way: This is a classical idea that, in a system of reaction-diffusion 
equations with a constant solution having the DDI property, one expects stable patterns 
to appear around that constant steady state. For the initial-boundary value problem 
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for reaction-diffusion-ODE system with a single diffusion equation (1.1)-(1.3), station- 
ary solutions can be constructed in the case of several interesting models (see Section 
3). However, the same mechanism which destabilizes constant solutions of such models, 
destabilizes also non-constant solutions. 

We prove Theorem 2.1 in Section 4 by applying ideas invented for analysis of the Euler 
equation and other fluid dynamics models. In that approach, it suffices to show that the 
spectrum of the linearized operator 

u\ = ( \ + ( M u > v ) Uu,v) \(u 

v J ~\Av J + y g u (U,V) g v (U,V) ) \ v 

with the Neumann boundary condition d v v = 0, has so-called spectral gap, namely, there 
exists a subset of the spectrum cr(£), which has a positive real part, separated from zero. 
Here, we prove that cr(£) C C consists of the interval {f u (U(x),V(x)) : x G Q} and of 
isolated eigenvalues of C, see Section 4 and, in particular, Fig. 4.1 for more details. 

More precise description of o~(£) is possible in the case of our model examples gath- 
ered in Section 3, where the autocatalysis condition (2.7) is satisfied, and moreover, the 
quantity f u (U(x), V(x)) is independent of x for every regular stationary solution. 

Corollary 2.4. Let (U, V) be a regular solution of problem (2.1)-(2.3). Assume that there 
exists a constant Ao > such that the following autocatalysis condition is satisfied: 

(2.9) < A = fu(U{x), V(x)) for all x G H. 

Suppose also the following "compensation condition" at the stationary solution (U,V) 

(2.10) g v (U(x),V(x))f v (U(x),V(x)) < for all xGfi. 

Then, the spectrum o~(C) contains the number X Q , which is an element of the continuous 
spectrum of £ and a sequence of real eigenvalues {X n }^ =1 of C converging towards Ao- 

In Section 3, we recall models from mathematical biology where both conditions: auto- 
catalysis (2.7) and compensation (2.10) are fulfilled, which can be checked in a surprisingly 
simple way. 

Remark 2.5. If one replaces the compensation condition (2.10) by the following one 

(2.11) g u (U(x),V(x))f v (U(x),V(x)) > for all x G U, 

then the statement of Theorem 2.4 remains true, what can be proved by an obvious 
modification of the proof given in Section 4. However, below in Corollaries 2.14 and 
2.15, we show that the compensation condition (2.10) has to be satisfied for reaction- 
diffusion-ODE problems which have constant solutions with the DDI property as well as 
non-constant stationary solutions. 

The initial-boundary value problem (1.1)-(1.4) may also have non-regular steady states 
in the case when the equation f(U,V) = is not uniquely solvable. Choosing different 
branches of solutions of the equation f(U(x),V(x)) = 0, we obtain the relation U(x) = 
k{V(x)) with a discontinuous, piecewise C 1 -function k. Here, we recall that a couple 
(U,V) G L°°(Q) x W 1 ' 2 (fl) is a weak solution of problem (2.1)-(2.3) if the equation 
f(U(x), V(x)) = is satisfied for almost all x G f2 and if 



VV(x) ■ V(p(x) dx + I g(U(x),V(x))(f(x) dx = 
for all test functions (p G W l,2 {Vt). 
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In this work, we do not prove the existence of such discontinuous solutions and we refer 
the reader to classical works [2, 19, 30] as well as to our recent paper [14, Thm. 2.9] for 
explanations how to construct such solutions to one dimensional problems using the phase 
portrait analysis. Our goal is to formulate a counterpart of the autocatalysis condition 
(2.7) which leads also to instability of weak (including discontinuous) stationary solutions. 

The following two corollaries can be proved in the same way as Theorem 2.1. 

Corollary 2.6. Assume that (U,V) is a weak bounded solution of problem (2.1)-(2.3) 

satisfying the following counterpart of the autocatalysis condition 

(2.12) Range f u (U,V) = {f u (U(x) J V(x)) :iGfi}c [A ,A ] 

for some constants < Ao < Ao < oo. Suppose, moreover, that there exists xq G f2 such 
that f u (U,V) is continuous in a neighborhood of xq. Then, (U,V) is an unstable solution 
the initial-boundary value problem (1.1)-(1.4). 

In fact, Theorem 2.1 is a particular case of Corollary 2.6, however, we prefer to state 
those two results separately to emphasize the role of the autocatalysis condition (2.7) in 
the proof of instability of stationary solutions. 

The following corollary covers the case when the generalized autocatalysis condition 

(2.12) fails to hold. 

Corollary 2.7 (Instability of weak solutions). Assume that the nonlinear term in equation 
(1.1) satisfies f(0,v) = for all v G M. Suppose that (U,V) is a weak bounded solution 
of problem (2.1)-(2.3) with the following property: there exist constants < Ao < A < oo 
such that 

(2.13) A < f u (U{x),V(x)) < A for all x E tt, where U(x)^0. 

Suppose, moreover, that there exists x G f2 such that U(x ) ^ and that the functions 
U — U(x) as well as f u (U, V) are continuous in a neighborhood of x . Then, (U, V) is an 
unstable solution the initial-boundary value problem (1.1)-(1.4). 

Remark 2.8. A typical nonlinearity satisfying the assumptions of Corollary 2.7 has the 
form f(u,v) = r(u,v)u and appears in the models where the unknown variable u evolves 
according to the Malthusian law with a growth rate r depending on u and other variables 
of the model. 

We defer the proofs of Theorems 2.1 and of Corollaries 2.4, 2.6, and 2.7 to Section 4. 

2.2. Turing-type mechanism of pattern formation. Let us recall that, in initial- 
boundary value problems for react ion- diffusion equations modeling biological pattern for- 
mation, 

(2.14) u t = eAu + f(u,v), v t = DAv + g(u, v), 

the nonlinearities are typically satisfying the inequalities f u > and g v < 0. Then, the 
DDI phenomenon requires the condition f u g v — f v g u > 0, and therefore, /„ and g u satisfy 
also the compensation condition f v g u < 0. Consequently, the signs of /„ and g u should 
be opposite and there are only two possibilities 

fu fv 



9u 9v 

The former on the right-hand side above is called an activator-inhibitor system, and the 
latter is a resource- consumer system. We refer the reader to the Murray books [21, 22] 
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and to the review article [33] as well as to references therein for more information on DDI 
in the case of two component react ion- diffusion systems and to the paper [31] in the case 
of several component systems. A detailed discussion of the DDI phenomena in the case 
of three component systems with some diffusion coefficients equal to zero can be found in 
the recent work [1]. 

Our next goal is to show that the autocatalysis inequality f u > and the compensation 
inequality f v g u < at a stationary solution (cf. the assumptions of Corollary 2.4) appear 
indeed in system (2.14) with e = exhibiting DDI if there exist spatially heterogeneous 
stationary solutions. 

Since it is natural to expect the existence of such steady states in models based on the 
Turing idea, let us first state a simple but fundamental property of stationary solutions 
to (1.1)-(1.4). 

Proposition 2.9. Assume that (U, V) is a non-constant regular solution of stationary 
problem (2.1)-(2.3). Then, there exists xq G fl, such that vector (u,v) = (U(xq),V(xo)^ 
is a constant solution of problem (2.1) -(2.3). 

The proof of Proposition 2.9 is short. It suffices to integrate equation (2.2) over Q and 
to use the Neumann boundary condition (2.3) to obtain J n g(U(x), V(x)) dx = 0. Hence, 
there exists xq £ f2 such that g(U(xo),V(xo)) = 0, because U and V are continuous. 
Thus, by equation (2.1), we also have f(U(xo),V(xo)) = 0. 

Next, we formulate our standing assumption on constant solutions. In this work, we 
consider only non- degenerate constant stationary solution (u, v) of the react ion- diffusion- 
ODE system (1. l)-( 1.3) , vectors i.e. (u,v) satisfying f(u,v) = and g(u,v) = 0. It 
means 

(2.15) f u (u,v)+g v (u,v)^0, det f /«(«»«)^ 0) and f u (u,v) ? 0. 

Notice that the first two conditions in (2.15) allow us to study the asymptotic stability of 
(u, v) as a solution to the corresponding kinetic system of ordinary differential equations 

du dv 

(2.16) — = f{u,v), —=g(u,v), 

by analyzing eigenvalues of the corresponding linearization matrix, see Remark 5.2 below 
for more details. On the other hand, by last condition in (2.15), the equation f(U, V) — 
can be uniquely solved with respect to U in the neighborhood of (u, v). 

In the case described by Proposition 2.9, we say that a non-constant solution (U, V) 
touches a constant solution (u,v). Now, we prove an important property of the con- 
stant solutions which are touched by non-constant solutions, and, for simplicity of the 
exposition, we limit our study to regular solutions. 

Theorem 2.10. Let (U(x),V(x)j be a regular stationary solution of problem (1.1)— (1.3) 
and assume that all constant stationary solutions which are touched by (U,V), are non- 
degenerate (cf. (2.15)). Then, at least at one of those constant solutions, denoted here by 
(u,v), the following inequality holds 

(2.17) TT^detf f fA\ f fA\ ) > 0. 

f u (u,v) V 9u{u,v) g v {u,v) J 

The proof of Theorem 2.10 is based on properties of solutions to the elliptic Neumann 
problem (2.4)-(2.5) (see Theorem 5.1, below) and we present it in Section 5. 
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Now, let us derive properties of a constant solution (u, v) which satisfies inequality 
(2.17). 

Remark 2.11. One should notice in the following corollaries that we do not impose the 
autocatalysis condition (2.7) to prove the instability of a constant steady state. Here, this 
condition is a consequence of the fact that the constant solution (u, v) is touched by a 
non-constant solution and because (u,v) satisfies inequality (2.17). 

Corollary 2.12. Denote by (u,v) a non-degenerate constant solution of problem (1.1)- 
(1.4), which satisfies inequality (2.17). Then, we have the following alternative: 

• either (u,v) is the constant solution of the reaction- diffusion- ODE system (1.1)- 
(1.4) with the DDI property, 

• or (u,v) is an unstable solution of the kinetic system (2.16). 

In the example discussed in Subsection 3.2, we show that both cases from Corollary 2.12 
can indeed occur in one model, depending on its parameters. 

Now, we have to emphasize the following necessary condition for the DDI property, 
which is obtained immediately from Corollary 2.12. 

Corollary 2.13. If problem (1.1)-(1.4) has only one constant solution (u,v) which is 
asymptotically stable as a solution to kinetic system (2.16) and if there exists a noncon- 
stant regular stationary solution which touches (u,v), than the constant steady state (u,v) 
is an unstable solutions of (1.1)-(1.4) (thus, (u,v) has the DDI property) . 

Next, we discuss the compensation condition (2.10) to link Theorem 2.10 with the the 
assumptions of Corollary 2.4. 

Corollary 2.14. Let (u,v) be a constant solution of problem (1.1)-(1.4) satisfying in- 
equality (2.17). If(u,v) is an asymptotically stable solution to kinetic system (2.16), then 
it satisfies the compensation inequality f v (u,v)g u (u,v) < 0. 

Thus, we obtain the following local version of the compensation assumption (2.10) for 
a regular non-constant stationary solution. 

Corollary 2.15. Let (U, V) be a non-constant regular stationary solution of problem 
(1. 1) -(1.4) . Denote by (u,v) a constant solution which is touched by (U,V) and which 
satisfies inequality (2.17). Assume that (u, v) is an asymptotically stable solution of kinetic 
system (2.16). Then, there exists x G Q such that f v (U(xo),V(x ))g u (U(x ),V(x )) < 0. 

In other words, every non-constant solution (U(x), V(x)) which touches a constant 
solution with the DDI property has to satisfy compensation property (2.10) at least at 
the point xq and, by continuity, in a neighborhood of xq. Notice, however, that we require 
this condition in Corollary 2.4 to be satisfied for all igO. 

3. Model examples 

In this section, our results are illustrated by applying them to nonlinearities from math- 
ematical biology. 

3.1. Gray-Scott model. Instability phenomena described in Section 2 appear, for exam- 
ple, in the following initial-boundary value problem with nonlinearities as in the celebrated 
Gray-Scott system describing pattern formation in chemical reactions [7] 

(3.1) Ut = -(B + k)u + u 2 v for x E H, t > 0, 

(3.2) v t = Av- u 2 v + B(l -v) for iGfl, t > 0, 
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where B and k are positive constants, with the zero-flux boundary condition for v and 
with nonnegative initial conditions. Here, every regular positive stationary solution (U, V) 
of the Neumann boundary-initial value problem for equations (3.1)-(3.2) has to satisfy 
the relation U = (B + k)/V, hence, 

(B + k) 2 

(3.3) AV-BV- y ' +B = for x eQ, 

(3.4) d u V = for xedQ. 

All continuous positive solutions of such boundary value problem in one dimensional case 
have been constructed in our recent paper [14, Sec. 5]. A construction of discontinuous 
stationary solutions to the react ion-diffusion-ODE problem for (3.1)-(3.2) can be also 
found in [14, Thm. 2.9]. 

Our instability results, formulated in Section 2.1, implicate that all stationary solutions 
(constant, regular as well as discontinuous) of the reaction-diffusion-ODE problem for 
equations (3.1)-(3.2) are unstable under heterogeneous perturbations. For the proof, it 
suffices to notice that the autocatalysis assumption (2.7) holds true because, for U = 
(B + k)/V, the function f u (U(x), V(x)) is independent of x and satisfies 

A = f u (U(x),V(x)) = -(B + k) +2U(x)V(x) = B + k > for all iGfi. 

This is the case of nonlinearities to which one can apply Corollary 2.1, because the com- 
pensation assumption (2.10) at a positive solution (U, V) is also valid due to the following 
calculation f v (U, V) ■ g u (U, V) = —2U 3 V < for all positive stationary solutions. 
Now, it suffices to apply either Theorem 2.1 or Corollary 2.4. 



3.2. Activator-inhibitor system. Now, we consider the following activator-inhibitor 
system with no diffusion of activator 

for x E U, t > 0, 
for x E Q, t > 0, 

supplemented with positive initial data and with the zero-flux boundary condition for 
v = v(x, t). We assume that r > and exponents satisfy p > 1, q, r > and s > 0. Such 
nonlinearities were proposed by Gierer and Meinhardt [6] to describe a system consisting 
of a slowly diffusing activator and a rapidly diffusing inhibitor. In their model, a high 
concentration of activator is supposed to induce morphogenetic changes in a biological 
tissue. 

First, we observe that every positive regular stationary solution satisfies the relation 
U p ~ l = V q , thus, the function V = V(x) is a solution of the boundary value problem 

(3.7) AV -V + V Q = for xGfi, 

(3.8) d u V = Q for x G <9f2, 

where Q = qr/(p— 1) — s. This problem with Q > 1 was considered in a series of papers, 
e.g. [35, 11, 25, 26]. In particular, when 1 < Q < g±| for N > 3, and 1 < Q < oo for 
N = 1, 2, it was shown that it has a positive solution V for a sufficiently large domain Q. 
We refer the reader to the review article [23] for more comments and references on this 
problem. 



(3.5) 
(3.6) 



Ut 



-u + 



Av 



v q 



U 
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Here, the autocatalysis condition (2.7) holds true for regular stationary solutions of 
system (3.5)-(3.6), because 

7TP-1 

/ u (?/,y) = -l+p— = -l+p>0 for t/P- 1 ^^. 

Notice that Ao = f u (U(x), V(x)) is independent of x, as required in Corollary 2.4. Next, 
we check directly the compensation assumption (2.10) 

/ TTP \ f TJ r ~ 1 \ TJP+r-l _ 

f v {U 1 V)-g u {U 1 V) = {-q^)[r—}=^q 1 ^<Q for all x G Q, 

for every positive stationary solution (U, V). As a consequence either of Theorem 2.1 or of 
Corollary 2.4, we obtain that all positive regular stationary solutions to system (3.5)-(3.6) 
with the Neumann boundary condition for v are unstable. 

By Proposition 2.9, every regular stationary solution of the Neumann boundary-initial 
value problem for system (3.5)-(3.6) has to touch the only constant positive solution 
(u,v) = (1,1) of system (3.5)-(3.6). In the case of the kinetic system corresponding to 
(3.5)-(3.6), it is shown in [24] that this constant solution is 

(1) asymptotically stable if the parameter r > is small; 

(2) unstable if the parameter r > is large. 

Thus, by Corollary 2.12, the constant solution (u, v) = (1, 1) of the reaction-diffusion-ODE 
system (3.5)-(3.6) with small r > 0, has the DDI property. 

3.3. Model of early carcinogenesis. The main motivation for the research reported in 
this work comes from our studies of the reaction-diffusion system of three ordinary/partial 
differential equations modeling the diffusion-regulated growth of a cell population of the 
following form 

(3.9) ut=[-^- dc]u for x G H, t > 0, 

\ u + v / 

(3.10) v t = —dbV + u 2 w — dv for iGfl, t > 0, 

(3.11) w t = DAw — dgiu — u 2 w + dv + k for x G f2, t > 0, 

supplemented with zero-flux boundary conditions for the function w and with nonnegative 
initial conditions, [14]. Here, the letters a,d c ,db,d g ,d, D, k denote positive constants. 

In the present paper, we focus on a reduced two-equation version of model (3.9)-(3.11), 
obtained by quasi-steady state approximation of the dynamics of v. Applying the quasi- 
steady state approximation in equation (3.10) (namely, assuming that v t = 0), we obtain 
the relation v = u 2 w/(db + d), which after substituting into remaining equations (3.9) and 

(3.11) yields the following initial-boundary value problem for the reaction-diffusion-ODE 
system 

(3.12) u t — ( d c )u for x G H, t > 0, 



1 + uw 



d 



(3.13) w t = DAw — d„w ; -u w + k for x G f2, t > 0. 

db + d 

Here, the autocatalysis assumption (2.9) and the compensation assumption (2.10) are 
satisfied by simple calculations, which are analogous to those in previous examples. As 
a consequence, all nonnegative stationary solutions of system (3.12)-(3.13) are unstable 
due to Theorem 2.1 and Corollary 2.7. This gives a counterpart of our results on three 
equation model (3.9)-(3.11) proved in [14]. 
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Numerical simulations suggest that the two-equation model exhibits qualitatively the 
same dynamics as (3.9)-(3.11). A rigorous derivation of the two equation model (3.12)- 
(3.13) from the model (3.9)-(3.11) as well as other properties of solutions to (3.12)-(3.13) 
are presented in Appendix A of this work. Our detailed studies of stability properties 
of space homogeneous solutions of the two equation model (3.12)-(3.13) are reported 
in Appendix B. In particular, by Corollary 2.12, constant steady states of (3.12)-(3.13) 
are either unstable solutions of the corresponding kinetic system or they have the DDI 
property. 



4.1. Existence of solutions. We begin our study of properties of solutions to the initial- 
boundary value problem (1.1)— (1.4) by recalling the results on local- in-time existence and 
uniqueness of solutions for all bounded initial conditions. 

Theorem 4.1 (Local-in-time solution). Assume that u ,v £ L°°(Q). Then, there exists 
T = T(||mo||oo, Ikolloo) > such that the initial-boundary value problem (1.1)— (1.4) has a 
unique local-in-time mild solution u, v £ L°°([0,T], L°°(f2)). 

We recall that a mild solution of problem (1.1)— (1.4) is a couple of measurable functions 
u, v : [0, T] x Q i— > R satisfying the following system of integral equations 



where e is the semigroup of linear operators generated by Laplacian with the Neumann 
boundary condition. Since our nonlinearities / = f(u, v) and g = g(u, v) are locally 
Lipschitz continuous, to construct a local-in-time unique solution of system (4.1)-(4.2), it 
suffices to apply the Banach fixed point theorem. Details of such a reasoning and the proof 
of Theorem 4.1 in a case of much more general systems of reaction-diffusion equations can 
be found in [29, Thm. 1, p. Ill], see also our recent work [14, Ch. 3] for a construction 
of nonnegative solutions of particular react ion- diffusion- ODE problems. 

Remark 4.2. If uq and v$ are more regular, i.e. if for some a £ (0, 1) we have Uq £ C a (Q), 
vq £ C 2+a (Q) and d u v = on dQ, then the mild solution of problem (1.1)— (1.4) is smooth 
and satisfies u £ C 1 ' a ([0,T] x H) and v £ C 1+a/2 ' 2+a ([0,T] x H). We refer the reader to 
[29, Thm. 1, p. 112] as well as to [5] for studies of general react ion-diffusion-ODE systems 
in the Holder spaces. 

4.2. General instability result. The next goal in this section is to show that, under ap- 
propriate conditions, linear instability of stationary solutions of a reaction-diffusion-ODE 
problem implies their nonlinear instability. Such a theorem is well-known for ordinary 
differential equations. Furthermore, in the case of react ion- diffusion equations, one can 
apply the abstract result from the book by Henry [8, Thm. 5.1.3]. However, in the 
case of reaction-diffusion-ODE problems, we often encounter discontinuous stationary so- 
lutions and the corresponding linearization operator may have a non-empty continuous 
spectrum. Hence, checking the assumptions of the general result from [8] does not seem 
to be straightforward. 

Therefore, we propose here a different approach which was developed to study stability 
of stationary solutions of fluid dynamic equations [10, 4], using the, so-called, linearization 



4. Instability of stationary solutions 



(4.2) 



(4.1) 
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principle. In that setting, only the existence of a spectral gap of a linearization operator is 
required to show the instability of steady states. This classical method has been recently 
recalled by Mulone and Solonnikov [20] to show the instability of regular stationary solu- 
tions to certain react ion-diffusion-ODE problems, however, assumptions imposed in [20] 
are not satisfied in our case. 

The crucial idea underlying this approach is to use two Banach spaces: a "large" space 
Z where the spectrum of a linearization operator is studied and a "small" space X C Z 
where an existence of solutions can be proved. 

Remark 4.3. To study instability properties of solutions to the reaction-diffusion-ODE 
problem (1.1)-(1.4), we use X = x and Z = L 2 (fi) x L 2 (Q) for a bounded 

domain Q C R , supplemented with the usual norms. 

Let us consider a general evolution equation 

(4.3) w t = £w + AT{w), w{0) = w 

where £ is a generator of a Co-semigroup of linear operators {e tc }t>o on a Banach space 
Z and M is a nonlinear operator such that A/"(0) = 0. 

Definition 4.4. Let (X, Z) be a pair of Banach spaces such that X C Z with a dense 
and continuous embedding. A solution w = of the Cauchy problem (4.3) is called 
(X, Z)-nonlinearly stable if for every e > 0, there exists 5 > so that if w(0) G X and 
\\w(0)\\ z < 5, then 

(1) there exists a global in time solution to (4.3) such that w G C([0, oo); X); 

(2) \\w(t)\\ z < e for allt G [0,oo). 

An equilibrium w = that is not stable (in the above sense) is called Lyapunov unstable. 

In this work, we drop the reference to the pair (X, Z). Let us also note that, under this 
definition of stability, a loss of the existence of a solution to (4.3) is a particular case of 
instability. 

Now, we recall a result linking the existence of the so-called spectral gap to the nonlinear 
instability of a trivial solution to problem (4.3). 

Theorem 4.5. We impose the two following assumptions. 

• The semigroup of linear operators {e tc }t>o on Z satisfies "the spectral gap condi- 
tion", namely, we suppose that for every t > the spectrum a of the linear operator 
e tc can be decomposed as follows: a = cr(e <£ ) = a_ U a + with a + ^ 0, where 

a_c{zeC : e Kt < \z\ < e#} and a + C {z G C : e Mt < \z\ < e At } 

and 

-oo < k < [i < M < A < oo for some M > 0. 

• The nonlinear term M satisfies the inequality 

(4.4) < Co||io||x||HU for all wEX satisfying \\w\\x < p 
for some constants Cq > and p > 0. 

Then, the trivial solution wq = of the Cauchy problem (4.3) is nonlinearly unstable. 

The proof of this theorem can be found in the work by Friedlander et al. [4, Thm. 2.1]. 
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Remark 4.6. The operator £ considered in this work satisfies the "spectral mapping 
theorem": a(e tc ) \ {0} = e ta( - c \ see Lemma 4.8, below. Thus, due to the relation \e z \ = 
e Rez for every z G C, the spectral gap condition required in Theorem 4.5 holds true if for 
every A G cr(£), either Re A G (k, fi) or Re A G (M, A). 

Remark 4.7. The authors of the reference [4, Thm. 2.1] formulated their instability result 
under the spectral gap condition for a group of linear operators {e tc } t ^R and in the case 
of a finite constant k (caution: in [4], the letter A is used instead of k). However, the 
proof of [4, Thm. 2.1] holds true (with a minor and obvious modification) in the case of a 
semigroup {e tc }t>o as well as k = — oo is allowed, as stated in Theorem 4.5. This extension 
is important to deal with the operator £ introduced in Lemma 4.8 below, which generates 
a semigroup of linear operators, only, and which may have an unbounded sequence of 
eigenvalues. We refer the reader also to the paper [32] for another instability result under 
less restrictive assumption on the spectrum of an operator £, but for a smaller class of 
nonlinearities. 

4.3. Linearization of reaction-diffusion-ODE problems. Let (U, V) be a stationary 
solution of problem (1.1)-(1.4). Substituting 

u = U + u and v = V + v 

into (1.1)-(1.2) we obtain the initial-boundary value problem for the perturbation (u,v) 
of the form (4.3): 

<«> £(?H(IM*' 

with the Neumann boundary condition, d v v = 0, where the linear operator £ and the 
nonlinearity M are defined and studied in the following two lemmas. 

Lemma 4.8. Let (U,V) be a bounded (not necessarily regular) stationary solution of 
problem (1.1)-(1.4). We consider the following linear system 

U * (ut\_ r (u\ = (u\( f u (U,V) f v (U,V) \fu 

[ } \v t J -^{v J ~{Av J + { g u (U,V) g v (U,V) J\v 

with the Neumann boundary condition d v v = 0. Then, the operator £ with the domain 
D(£) = L 2 (Q) x W 2 ' 2 (Q) generates an analytic semigroup {e tc } t >o of linear operators on 
L 2 (Q) x L 2 (Vt), which satisfies "the spectral mapping theorem": 

(4.7) a(e tc ) \ {0} = e ta(c) for every t>0. 

Proof. Notice that £ is a bounded perturbation of the operator 

T* W o 
v Au 



£r 



with the Neumann boundary condition for v and with the domain D(£q) = L 2 (Q) x 
W 2,2 (Q), which generates an analytic semigroup on L 2 (Q) x L 2 (Q). Thus, it is well-known 
(see e.g. [3, Ch. III. 1.3]) that the same property holds true for £. 

"The spectral mapping theorem" for the semigroup {e tC } t >o expressed by equality (4.7) 
holds true if the semigroup is e.g. eventually norm-continuous (see [3, Ch. IV. 3. 10]). Since 
every analytic semigroup of linear operators is eventually norm-continuous, we obtain 
immediately relation (4.7) (cf. [3, Ch. IV, Corollary 3.12]). □ 

Next, we show that the nonlinearity in equation (4.5) satisfies the assumption (4.4) 
from Theorem 4.5. 
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Lemma 4.9. Let (U,V) be a bounded (not necessarily regular) stationary solution of 
problem (1.1)-(1.4). Then, the nonlinear operator 

u\ = ff(U + u,V + v)-f(U,V)\_f f u (U,V) f v (U,V) \fu 
v J ~ I g(U + u,V + v)-g(U,V) J \g u (U,V) g v (U,V) J \v 



satisfies 



W(u,v)\\ L2xL2 < C(p, \\U\\ L ~, \\V\\ L ~)(\\u\\ L ~ + ||5T|U~)(||Ti|| £a + lluHia) 
= C(p, \\U\\ L oo, \\V\\ L ^)\\(u,v)\\ L oo xL oo\\(u,v)\\ L 2 xL 2 

for all u, v £ L°° such that \\u\\l°° < p and \\v\\l°° < p, where p > is an arbitrary 
constant. 

Proof. This is nothing but the Taylor formula applied to the C 2 -nonlinearities f = f(u, v) 
and g = g{u,v) in problem (1.1)-(1.2). □ 

4.4. Instability of steady states. Now, we are in a position to apply the described 
above general theory to show the instability of stationary solutions to react ion- diffusion- 
ODE problems. 

Proof of Theorem 2.1. Let (U(x), V(x)) be a regular stationary solution to problem (1.1)- 
(1.4). To show its instability, in view of Theorem 4.5 combined with Lemmas 4.8 and 
4.9, it suffices to study the spectrum a(£) of the linear operator £ defined formally by 
formula (4.6) with the domain D(£) = L 2 (Q) x W 2 ' 2 (Q). 
Let us define the constants 

(4.8) A = mlf u (U(x),V(x)) > and A = sup /„(£%), V(xj) > 0, 

where the positivity of Ao is a consequence of the autocatalysis condition (2.7). We are 
going to prove that cr(£) C C consists of all numbers from the interval [Ao,Ao] and of a 
set of (possibly complex) eigenvalues of (£,D(£)) which are isolated points of C. 
Part I: Interval [Ao, Aq]. First, we show that for each A £ [Ao, Aq] the operator 

£ - XI : L 2 (n) x W 2 > 2 (n) -> L 2 (Sl) x L 2 {Sl) 

defined by formula 

{£- XI)((p,ip) = ((f u - \)<p + f v il>, + g u <f+ (g v - A)V>), 

where f u = f u (U(x), V(x)), etc., cannot have a bounded inverse. Suppose, a contrario, 
that (£ — A/) -1 exists and is bounded. Then, for a constant K = \\(£ — A/)" 1 )!, we have 

IKv ? ; V ; )IU 2 (n)xw2.2(f7) < K\\(£ - \I)(<p,il>)\\L*(n)xL2{n) 

for all (if, ip) £ L 2 (Vl) x W 2,2 (VL) or, equivalently, using the usual norms in L 2 {VL) x W 2 ' 2 (Q) 
and in L 2 (Q) x L 2 (Q): 

IklU 2 ^) + IMlHW(n) 

(4.9) 

< K(\\(fu - A)v? + M\\v[a) + || + 9u^ + (g v ~ A)^|| L 2(fi))- 

A contradiction will be obtained by showing that inequality (4.9) cannot be true for all 
M) £ L 2 (tt) x W 2 ' 2 (tt). 

To prove this claim, first we observe that, for each A £ [Ao,Ao], there exists xq £ Q 
such that f u (U(xo), V(xq) \ — A = 0. Hence, for every e > there is a ball B £ C Vt such 
that ||/ u - A|| L oo (Be) < e. 
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Next, for arbitrary xjj G C^°(Q) such that supp?/> C B £ , we choose tp G L 2 (Q) such 
that supp tp C B £ and in such a way that Aip + g u p + (g v — \)ip = £, where the function 
( G L 2 (fl) satisfies ||C||,L 2 (n) < £|MU 2 (fi)- Let us explain that such a choice of <p, ( G L 2 (Q) 
is always possible. We cut g u at the level e in the following way 



9l 



nt(TT(<r\ V(t\\ = S 9u(U(x),V(x)) if \gJU(x),V(x))\ >s, 
g u {U{x),v[x)) _ i £ if |^ tt (i7(a;),y(ar))|<e. 



Since, A^ + + (g v - A)^ = A^ + + (9v ~ A)^> + (g u - g e u )p, we may set 
(A^ + fa-A)^) ^ f= _ (9o _ s>£L2(!i) 

with - ^Hloo(q) < e. 

Now, noting that supp p C we obtain the inequality 

||(/« - A)^|| L 2 (n) < - A) 11^(^)11^1^2^) < e\\p\\ L 2 {n) . 

Thus, substituting functions p, ip, and ( into inequality (4.9), we obtain the estimate 

IMU 2 (n) + ||^||w/ 2 . 2 (n) 

(4.10) < K{\\(f u - X)\\ L ~ m \\p>\\ L2m + ||/WlU 2 (n) + KWmn)) 
< K(2e\\ip\\ L 2 {n) + \\f v \\L°°(n)\\tp\\L2(n))- 

Hence, choosing e > in such a way that 2Kb < 1 and compensating the term 2Kb \\<p\\l 2 (o.) 
on the right-hand side of inequality (4.10) by its counterpart on the left-hand side, we 
obtain the estimate llV'lliy 2 ' 2 ^) < -^||/'y|U o °(f2)||'0l|i 2 (n)) which, obviously, cannot be true 
for all ip G C^°(Q) such that supp^ C B £ . We have completed the proof that each 
A G [A , A ] belongs to cr(C). 

Part II: Eigenvalues. In the next step, we show that the remainder of the spectrum of 
(£,£>(£)) consists of a discrete set of eigenvalues {A„}^ =1 C C \ [A ,A ], analyzing the 
corresponding resolvent equations 

(4.11) (f u -\)cp + f v iP = F in H 

(4.12) Aip + g u <p + (g v - X)i[) — G in Q 

(4.13) d u ip = on dtt, 

with arbitrary F,G G L 2 (Q). Here, one should notice that for every A G C \ [Ao, Ao], 
one can solve equation (4.11) with respect to <p. Thus, after substituting the resulting 
expression <p — (F — f v ip)/(f u — A) G L 2 (Q) into (4.12), we obtain the boundary value 
problem 

(4.14) Aip + q(X)i(j — p(A) for x G fi, 

(4.15) d v ip = for x G dQ, 

where 

(4.16) q(X) = q(x, A) = --^K + g v - A and p(X) = p(x, A) = G 



fu — A f u — X 

For a fixed A G C \ [Ao,A ], by the Fredholm alternative, either the inhomogeneous 
problem (4.14)-(4.15) has a unique solution (so, A is not an element of &(£)) or else the 
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spectral gap 



\ 



A A /uo 
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Figure 4.1. The spectrum a(£) is marked by thick dots and by the inter- 
val [Ao, Ao] in the sector T,$^ - The spectral gap is represented by the strip 
{A G C : fi < Re A < M} without elements of er(£). 

homogeneous boundary value problem 

(4.17) Atp + q(\)tp = for x G Q, 

(4.18) <9^ = for xedtt, 

has a nontrivial solution ip. Hence, it suffices to consider those A G C \ [Ao, Ao], for which 
problem (4.17)-(4.18) has nontrivial solution. 

Now, we prove that the set cx(£) \ [Ao,A ] consists of isolated eigenvalues of £, only. 
First, we rewrite problem (4.17)-(4.18) in the form 

1> = G[-(q(\)+e)1>]=R{\)4, 

where the operator G = "(A — £I)~ ln supplemented with the Neumann boundary con- 
ditions is defined in the usual way. Here, I G ffi is a fixed number different from each 
eigenvalue of Laplacian with the Neumann boundary condition. 

Recall that, for each A G C \ [A , A ], the operator -R(A) : L 2 (Q) — > L 2 (Q) is compact as 
the superposition of the compact operator G and of the continuous multiplication operator 
with the function g(A) +£ G L°°(Q). Moreover, the mapping A i— > -R(A) from the open set 
C \ [Ao, A ] into the Banach space of linear compact operators is analytic, which can be 
easily seen using the explicit form of g(A) in (4.16). Thus, the set ct(£) \ [A , A ] consists 
of isolated points due to the analytic Fredholm theorem (see Theorem 4.11, below). Here, 
to exclude the case (a) in Theorem 4.11, we have to show that the operator / — R(X) is 
invertible for some A G C \ [A , A ]. This is, however, the consequence of the fact that 
the inhomogeneous boundary value problem (4.14)-(4.15) has a unique solution if A > 
is chosen so large that q(x, A) < 0. 

Part III: Spectral gap. By Lemma 4.8, there exists a number uq > such that the 
operator (£ — uoI,D(£)) generates a bounded analytic semigroup on L 2 (Q) x L 2 (Q), 



16 



ANNA MARCINIAK-CZOCHRA, GRZEGORZ KARCH, AND KANAKO SUZUKI 



hence, this is a sectorial operator, see [3, Ch. II, Thm. 4.6]. In particular, there exists 
5 G (0,7r/2] such that a(£) C £ 5jW0 = {A G C : |arg (A - u )\ > tt/2 + 5}, see 
Fig. 4.1. The part of the spectrum a(£) in the triangle fl {A G C : Re A > 0} 
consists of all numbers from the interval [Ao, Ao] with Ao > and of a discrete sequence of 
eigenvalues with accumulation points from the interval [Ao, Ao], only. Thus, we can easily 
find infinitely many < \x < M < Ao, for which the spectrum o~(£) can be decomposed 
as required in Theorem 4.5. Here, one should use the spectral mapping theorem, i.e. 
equality (4.7), and Remark 4.6. □ 

Remark 4.10. One should emphasize here that the instability of steady states from The- 
orem 2.1 is caused not by an eigenvalue with a positive real part, but rather by positive 
numbers from the interval [Ao, Ao] which are, in some cases, elements of continuous spec- 
trum of the operator (£,D(£)), see Corollary 2.4. 

The following general result on a family of compact operators was used in the proof of 
Theorem 2.1. We state it here for reader's convenience and its proof can be found in the 
Reed and Simon book [28, Thm. VI. 14]. 

Theorem 4.11 (analytic Fredholm theorem). Assume that H is a Hilbert space and 
denote by L(H) the Banach space of all bounded linear operators acting on H . For an 
open connected set D C C, let f : D — >• L(H) be an analytic operator-valued function such 
that f(z) is compact for each z G D. Then, either 

(a) (J — f(z))' 1 exists for no z G D, or 

(b) (I — /(z)) _1 exists for all z G D\S, where S is a discrete subset of D (i.e. a set 
which has no limit points in D). 

A more precise description of cr(£) is possible in the case of particular nonlinearities, 
as stated in Corollary 2.4. First, however, we recall properties of eigenvalues /i and the 
corresponding eigenfunctions if) of the following boundary value problem 



with a given potential q G L°°(Q). The following two lemmas contain results on the exis- 
tence of eigenvalues, their continuous dependence on a potential q, and their monotonicity 
with respect to q. 

Lemma 4.12. Assume that q G L°°(Q). Then, there exists a sequence {fi>k(q)}kLo of 
eigenvalues of problem (4.19)-(4.20) satisfying fii(q) < /^((Z) < fail) +oo. 

Lemma 4.13. There exists a constant C > such that for all q,q G L°°(Q) and for all 
n G NU{0}, we have \fi n (q) — A*n(5)| < C\\q — q\\L 2 (n)- Moreover, if q(x) < q(x) a.e., then 

Results stated in these lemmas are classical. Indeed, by a well-known reasoning, it 
suffices to convert problem (4.19)-(4.20) into an eigenvalue problem for a compact, self- 
adjoint, nonnegative operator on L 2 (Q) (as in Part II of the proof of Theorem 2.1). Hence, 
properties of eigenvalues, stated in these lemmas, result immediately from the abstract 
theory, see eg. [27, Cor. 5.6 and Thm. 5.7]. 

Proof of Corollary 2.4- As in the proof of Theorem 2.1, we analyze the resolvent equations 
corresponding to £ stated in (4.11)-(4.13) with f u = X . 

It follows from equations (4. 11)- (4. 13) that the operator £ — X I : D(£) — > L 2 {Vt) x 
L 2 (Q) is injective. Indeed, setting A = Aq = f u and F = in equation (4.11) we obtain 



(4.19) 
(4.20) 



Atp + fiqip = 
d^ = 



for 

for x G dQ, 
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that ip = almost everywhere, because f v y^0 due to the compensation condition (2.10). 
Thus, using equation (4.12) with ip = and G = 0, we obtain that ip = because g u ^ 
by (2.10). Hence, Ker(£ - A I) = {(0,0)}. However, (£ - A /) _1 is not well-defined on 
the whole space L 2 (Q) x L 2 (Q), because, by equation (4.11) with A = Ao, the relation 
f v ip = F cannot be true for every F G L 2 (Q) and some ip G W 2,2 (Q). Thus, the number 
Ao = fu belongs to the continuous spectrum of (£, D(C)). 

Next, if A 7^ /„, as it was explained in the proof of Theorem 2.1, it suffices to analyze 
nontrivial solutions of the boundary value problem in Q 

A^ + g(A)^ = with q (\) = -J^. + 9v -\ 
d u tp = on <9f2. 

Step 1. Notice that the function q = q(X, •) : M\{Ao} — > L 2 (Q) is continuous and satisfies 
q(X,x) — > +oo uniformly in x G Q if A / A , because of the autocatalysis assumption 
g u fv < (here, q(X,x) — > +oo also if g u f v > and if A \ Ao, see Remark 2.5). 

Step 2. By Lemma 4.12, for every A < Ao, there exists an increasing sequence 

/ii(g(-, A)) < fi 2 (q(-, A)) < ... < Hn{q{-, A)) < ... -> +oo 

of eigenvalues of problem (4.19)-(4.20) with q = q(x, A). Our goal is to show that there 
exists a sequence of real numbers A n < Ao such that /i n (g(-, A n )) = 1. Then, the corre- 
sponding eigenfunction of problem (4.19)-(4.20) will be a non-zero solution of (4.21). 

Step 3. The eigenvalue [i n (q(-, A)) is a continuous function of A < Ao for each n G N. 
Indeed, this is an immediate consequence of Lemma 4.13, because the mapping q(-, A) : 
(— oo, Ao) — > L 2 (Q) is continuous. 

Step 4- Let us show that fi n (X) — > when A — > Ao- By Step 1, we have 

uj(X) = inf q(x, A) — > oo when A < Ao and A — > Xq. 

Since cu(A) is a positive constant for each A < A sufficiently close to A , then the eigen- 
values /i of the problem 

+ noj(X)ip = in VI 
d u i/j = on dil. 

has a sequence of eigenvalues /i n (a;(A)) > satisfying /i n (w(A)) — >■ as uj(X) — > oo for 
each n G N. Since uj(X) < q(x,X), the comparison property for eigenvalues stated in 
Lemma 4.13 implies 

/i n {q(-, A)) < n n (u)(\)) -)• as A -> A . 

Conclusion. Since < /ii(g(-,A)) < ^(^(^A)) < ... < fi n (q(-,X)) < ... — > +oo and 
X n (q(-, A)) — > as A — > Ao, for every n G N, it follows from the continuous dependence of 
/i n (g(-, A)) on A (see Step 3) that there exists A„ — > Ao such that fi n (q(-, X n )) = 1, provided 
n is sufficiently large. 

Thus, o~(C) contains the number Ao = f u > and a sequence of real numbers A n /• Ao- 
Notice, however, that the spectrum cr(C) may consist also of other complex eigenvalues 
of the operator C, according to the reasoning from the proof of Theorem 2.1. □ 

Proof of Corollary 2.6. Here, it suffices to follow the proof of Theorem 2.1. First, we 
linearize our problem at a weak bounded stationary solution (U, V) and we notice that 
assumptions of Lemmas 4.8-4.9 are satisfied. Next, following the arguments in Part I of 
the proof of Theorem 2.1, we show that the number f u (U(xo), V(xq)) belongs to cx(£), 
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where f u (U(x),V(x)) is positive at xq and continuous in its neighborhood. Notice that 
we do not need to show that all numbers from Range f u (U, V) are in <j(£) to show the 
spectral gap condition required by Theorem 4.5. The reasoning from Part II of the proof 
of Theorem 2.1 can be copied here without any change because q(X,x) defined in (4.16) 
is a bounded function for every A e C \ [Ao, A ]. □ 

Proof of Corollary 2. 7. Here, the analysis is similar to the case of regular stationary so- 
lutions discussed in Theorem 2.1, hence, we only emphasize the most important steps. 

Let (£/, V) be a weak solution of problem (2.1)-(2.3) and denote by X C ft its null set, 
namely, a measurable set such that U(x) =0 for all x G X and U(x) ^ for all x G ft \X. 
For a null set X, we define the associate L 2 -space 

L 2 (Q) = {v G X 2 (ft) : v(x) = for x G X}, 

supplemented with the usual L 2 -scalar product, which is a Hilbert space as the closed 
subspace of L 2 (ft). In the same way, we define the subspace L^(Q) C L°°(ft) by the 
equality Lf(Q) = {v G L°°(ft) : v(x) = for x G X}. 

Obviously, when the measure of X equals zero, we have X|(ft) = L 2 (ft). Assumptions 
imposed in Corollary 2.7 imply that X is different from the whole interval. 

Now, observe that if Uq(x) = for some x G ft then by equations (1.1) with the 
nonlinearity f(u,v) = r(u,v)u we have u(x,t) = for all t > 0. Hence, the spaces 

(4.22) X x = L^(p) x L°°(ft) and Z x = Lf (ft) x L 2 (tt) 

are invariant for the flow generated by problem (1.1)-(1.4) (notice that there is no "X" in 
the second coordinates of Xx and Zx)- The crucial part of our analysis is based on the 
fact that, as long as we work in the space Xx and Zx, we can linearize problem (1.1)-(1.4) 
at the weak solution (U,V). Moreover, for each x G ft \ X, the corresponding linearized 
operator agrees with £ defined in Lemma 4.8. Hence, the analysis from the proof of 
Theorem 2.1 can be directly adapted to discontinuous steady states in the following way. 

We fix a weak stationary solution (Ux, V%) with a null set X C ft. The Frechet de- 
rivative of the nonlinear mapping J 7 : Zx — >■ Zx defined by the mappings (U, V) (-> 
(f(U, V),g(U, V)) at the point (U x , V x ) G Z x has the form 

mu* vx) ( ^ ) = Ax(x) 

where 

AM-f fu(Ux(x),V x (x)) f v (U x (x),V x (x)) \ 
] ~ V 9u(U x (x),V x (x)) g v (U x (x),V x (x)) J 

This results immediately from the definition of the Frechet derivative. 
Next, we study spectral properties of the linear operator 

in the Hilbert space Z x (see (4.22)) with the domain D(C X ) = Lf(ft) x W 2 > 2 {Vt). Here, 
the reasoning from Part I of the proof of Theorem 2.1 can be directly adapted with a 
modification as in the proof of Corollary 2.6. 

Finally, we study the discrete spectrum of C x in the same way as in Part I of the proof 
of Theorem 2.1 because the corresponding function q(X, x) is bounded for A G C \ [A , A ]. 

The proof of instability of the stationary solution (U x , V x ) is completed by Theorem 4.5 
and Lemmas 4.8-4.9. □ 
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5. Constant and non-constant stationary solutions 
First, we prove an auxiliary result which will imply Theorem 2.10. 

Lemma 5.1. Assume that V G C 2 (Q)C\C 1 (Q) is a non- constant solution to the boundary 
value problem 

(5.1) AV + h(V) = for x en and d v V = for x E dtt. 
Then, there exists xq G Q and ao G R such that 

(5.2) V^(xo) = ao, h(ao) = and h'(ao) > 0. 

Proof. First, as in the proof of Proposition 2.9, we integrate the equation in (5.1) and we 
use the Neumann boundary condition to obtain J Q h(V(x))dx = 0. Hence, there exists 
x G Q and a G R such that V(xq) = ao and h(ao) = 0. Now, we suppose that h'(ao) < 0, 
and consider two cases: x G D, and x G dQ, separately. 
Let xo G Q. Since h(ao) = 0, we have 

A(V - a ) + h(V) - h(a ) =0. 

Using the well-known formula 



f 1 d 

I -j-h(sV + (1 - s)a ) ds 
Jo ds 

(V-a )[ h'(sV+(l-s)a )ds, 



h{V) - h(ao) 



we obtain 

(5.3) A{V-a ) + r(x,a )(V -a ) = 0, 

where r(x, ao) = h'{sV(x) + (1 — s)ao) ds. Observe that r(-, ao) G C{Vt) and 

r(x ,a )= / h'(sV(x ) + (1 - s)a ) ds 



h'(sao + (1 — s)ao) ds 
= h'(a ) < 0. 

Hence, there exists an open neighbourhood U C Vt of Xo such that r(x, a ) < for all 
x G U. Suppose that r(x, a ) < for all iffl. Multiplying both sides of equation (5.3) 
by V(x) — a and integrating over Q, we obtain 

\V(V(x) — a )\ 2 dx + / r(x, a Q )(V(x) — a ) 2 dx = 0. 



This implies that V(x) = ao, which is a contradiction, because we assume that V = V(x) 
is a non-constant solution. Therefore, there exists x\ G dUC\Q such that r(x\,ao) = 0. It 
follows from equation (5.3) that AV(x\) = 0, and consequently, from equation (5.1) we 
have h{y(x\)) = 0. Hence, there exists ai G R such that V{x\) = ai and h(ai) = 0. Note 
that a\ 7^ a - Thus, if the equation h(V) = has only one root, the proof of Lemma 5.1 
is completed. 

Now, we consider two cases. 
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Case I: The equation h(V) = has no solution between a and a±. Thus, we define the 
function 

ip(s) = V(xo + s(x\ — xo)) for < s < 1, 

and, without loss of generality, we can assume that a$ < tjj(s) < ai for all < s < 1. 

Since h(ao) = and h'(a ) < 0, we have h(ao + 9 ) < for small 9 > 0. If we 
also suppose that h\a\) < 0, then, we can find small B\ > such that h(ai — 6{) > 0. 
Noting that ip(s) is continuous function, we see that there exist s*, s** G (0, 1) such that 
ip(s*) = a + 9o and tp(s**) = a± —6\. This implies that there exists s G (s*, s**) for which 
h(V(xo + s(xi — Xo))) = 0, and from the assumption for if}(s), 

a < V(x + s(x\ — Xo)) < di- 

This is a contradiction with the hypothesis that the equation h(V) = has no roots 
between a and a\. Hence, h'(ai) > 0. 

Case II: The equation h(V) = has a solution a m between and a±. It is clear that 
V(x) has to touch a m , too. Choosing a m the closest root of h(V) = to a , we repeat the 
argument from Case I to show that h'(a m ) > 0. 

Next, let Xo G dQ. Following the previous reasoning and using the hypothesis h'(ao) < 
0, we find a ball B C Q such that xq G dQ and r(x, ao) < for all x G B. Moreover, 
we can assume that either V(x) > a or V(x) < a for all x G B, because, if there exists 
x% G B such that V(x±) = ao, then we can apply the same argument as in the first part 
of this proof to obtain h'(ao) > 0. 

Let V(x) < a for all x G B, and we apply the Hopf boundary lemma to equation (5.3) 
in the ball B. If V is a non-constant solution satisfying V(x) — a < and V(x ) — a = 0, 
then necessarily dV(xo)/du > 0, which contradicts the Neumann boundary condition 
satisfied by V at xo G dQ. 

Now, we consider the case V(x) > ao for all x G B. Here, the function U(x) = 
— (V(x) — ao) satisfies the equation 

-AU + (-r(x,a ))U = mB 

where r(x,a ) < 0, U(x) < for all x G B and U(x ) = 0. Hence, the Hopf boundary 
lemma yields a contradiction. 

Thus, h'(a ) > and the proof is complete. □ 

Proof of Theorem 2.10. Let (U(x),V(x)) be a non-constant regular stationary solution 
of (1.1)— (1.3) which means that U = k(V) and f(k(V),V) = 0. Since each constant 
solution is non-degenerate, the equation f(U,V) = can be solved locally with respect 
to U, which implies that k = k(V) is a C 1 -function. Substituting U = k(V) into equation 
(2.2) and denoting h(V) = g[k(V), V), we obtain the following boundary value problem 
satisfied by V = V(x): 

(5.4) AV + h(V) =0 for x G Q, 

(5.5) d v V = for x G dtt. 

By Lemma 5.1, there exists a constant solution v of problem (5.4)-(5.5) and Xo G fl such 
that 

(5.6) V(x ) = v, h{v) = 0, and h'(v) > 0. 

On the other hand, differentiating the function h{v) = g(k{v),v) yields 

(5.7) h'{v) = k'(v)g u (k(v),v)+g v (k(v),v). 
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Moreover, we differentiate both sides of the equation f(k(v),v) = with respect to v to 
obtain k'(v)f u (k(v),v) + f v (k(v),v) = 0. Hence, 

Finally, choosing v = v and u = u = k{v) and substituting equation (5.8) into (5.7), we 
obtain 

, lf -s fv(k(v),v) ( , . _^ , . . — \ 

= 7 }- -x [/u(«, U) - uJfl'uCw, «)] 



det 



fu(u,v) fv(u,v) 



fu(u,v) \ g u {u,v) g v {u,v) 

By (5.6), it holds h'(v) > 0. Moreover, since (u,v) is non-degenerate, we obtain h'(v) > 0. 
This completes the proof of inequality (2.17). □ 

Remark 5.2. Before deriving consequences of inequality (2.17), let us recall properties of 
a constant solution (u, v) of the kinetic system of the ordinary differential equations 

(5.9) d i = fM ' Tt =9M - 
It is well-known that if 

(5.10) f u (u, v) + g v (u, v)<0 and det ( fu ^ V \ fv ^ V \ ) > 0, 



9u[u,v) g v {u,v) 



then the Jacobi matrix 

f u (u,v) f v (u,v) 
g u (u,v) g v (u,v) 



(5.11; 



has all eigenvalues with negative real parts, so, (u, v) is an asymptotically stable solution 
of system (5.9). On the other hand, if 

(5.12) either /.(«.»)+*(«,»)><> or det ( J K „_ 

then the linearization matrix (5.11) has an eigenvalue with a positive real part, and 
consequently, the couple (u,v) is an unstable solution of (5.9). 

Proof of Corollary 2.12. Let (u, v) be a constant solution of (1.1)-(1.4) for which inequal- 
ity (2.17) holds true. Using properties of constant solutions to the kinetic system recalled 
in Remark 5.2 we obtain immediately the following alternative. 

If both factors on the left-hand side of inequality (2.17) are positive, then, in particular, 
the autocatalysis condition f u (u,v) > is satisfied. Hence, the constant solution (u, v) is 
an unstable solution of the reaction-diflusion-ODE system (1.1)-(1.4) by Theorem 2.1. 

On the other hand, if both factors on the left-hand side of inequality (2.17) are negative, 
then the alternative in (5.12) is satisfied and the constant vector (u, v) is an unstable 
solution of the corresponding kinetic system (1.1). □ 
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Proof of Corollary 2.13. Since (u,v) is an asymptotically stable solution of the kinetic 
system (5.9), inequality (2.17) together with the second inequality in (5.10) imply im- 
mediately the autocatalysis condition f u (u,v) > 0. Hence, the instability of (u, v), as 
a solution of the react ion- diffusion- ODE problem (1.1)- (1.4), is a direct consequence of 
Theorem 2.1. □ 

Proof of Corollary 2.14- The vector (u,v) is an asymptotically stable solution of the ki- 
netic system (5.9), hence, inequalities (5.10) hold true (because we assume that (u,v) is 
nondegenerate, cf. (2.15)). Thus, calculating the determinant in the second inequality of 
(5.10) we obtain 

(5.13) f v (u, v)g u (u, v) < f u (u, v)g v (u, v). 

Next, by the first inequality in (5.10) and by Corollary 2.13 (see its proof), we have the 
following two inequalities 

f u (u,v) + g v (u,v) < and f u (u,v) > 0, 

which imply g v (u, v) < and, consequently, f u (u,v)g v (u,v) < 0. Using this estimate in 
inequality (5.13), we obtain the compensation property (2.10) at the constant steady state 
(u,v). □ 



Appendix A. Model reduction 

Initial-boundary value problems for reaction-diffusion-ODE systems arise in the mod- 
eling of the growth of a spatially-distributed cell population, where proliferation is con- 
trolled by endogenous or exogenous growth factors diffusing in the extracellular medium 
and binding to cell surface as proposed by Marciniak and Kimmel in the series of recent 
papers [15, 16, 17]. Here, we consider the following particular case of such systems, which 
was studied by us in [14] 

(A.l) d t u £ = (— ^ d c )u £ = f(u £ ,v £ ), 

(A.2) ed t v £ = -d b v £ + u £ 2 w £ - dv £ = g(u £ , v £ , w £ ), 

(A.3) d t w £ - DAw £ + d g w £ = -u £ 2 w £ + dv £ + k q = h(u £ , v £ , w £ ), 

on (0, oo) x Q, supplemented with zero-flux boundary conditions for w £ 

(A.4) d v w £ (x,t) = for all t > 0, x G dQ 

and with nonnegative, smooth and bounded initial data 

(A.5) u £ (x, 0) = u (x), v £ (x, 0) = v (x), w £ (x, 0) = w (x). 

In equations (A.1)-(A.3), the letters a, d c , db, d g , d, k , D denote positive constants. 

As it was shown in [14, Sec. 3], solutions of this problem are nonnegative and stay 
bounded on every interval [0, T\. Here, we study the behavior of these solutions as e — > 0. 

First, we notice that choosing e = in equation (A.2), we obtain the identity 



UNSTABLE PATTERNS IN REACTION-DIFFUSION-ODE SYSTEMS 23 

Substituting it to the two remaining equations (A.l), (A. 3) yields the system 

(A.7) u t =(-^- d c )u forxeH, t > 0, 

V 1 + uw J 

(A. 8) w t = DAw — d a w : — b ——u 2 w + k for x G f2, t > 0, 

d b + d 



(A.9) d v w(x, t) = for x G t > 0, 

(A. 10) w(x, 0) = Uq{x) and 0) = wq(x). 

Repeating our reasoning from [14] one can show that this new system has also a unique, 
global-in-time nonnegative, smooth solution for nonnegative and smooth initial data. 

In this part of Appendix, we show that solutions of the quasi- stationary system (A. 6)- 
(A.10) are good approximation of solutions of the original three-equation model (A.l)- 
(A.5). Our main result concerns uniform estimates for an approximation error for u and 
w on each finite time interval [0,T]. 

First, we show that solutions of e-problem (A.1)-(A.5) are uniformly bounded with 
respect to small e > 0, locally in time. 

Lemma A.l. For every T > there exists C(T) > such for all sufficiently small e > 
(e.g. e G (0, (d b + d)/(2d g ))) the solution of system (A.1)-(A.5) satisfies 

lk(*)IU~(n) < C(T), K(t)IU»(n) < C(T), K(t)IU~(n) < C(T) 
for allt G [0,T]. 

Proof. Since v £ /(u £ + v e ) < 1 for nonnegative solutions, equation (A.l) yields the inequal- 
ity, 

(A.ll) K(*)IU~(n) < M(t) = ||m |U-(o) 



(a—dc)t 



Hence, we have the estimate ||w e (i)||L°°(n) < C(T) = \\u \\L°°(n)^ a ~ dc ^ T for all t G [0, T\. 

Applying the comparison principle to the parabolic equation (A. 3) with the Neumann 
boundary condition we obtain the estimate 

(A.12) < w £ (x,t) < |K(0IU~ ( n) < C w (t), 

where the function C w = C w (t) satisfies the Cauchy problem 

^\ -fljp w + d g C w = rf||u E (*)||i°°(n) + K o 

C w (0) = \\w \\L°°(n), 

and is given by the formula 

(A.14) C w (t) = ^+ (\\w Q \\ Loom - e-^ + d j\- d ^\\v £ (r)\\ Lac{n) dT. 

Next, we use equation (A. 2) to obtain 

(A.15) ||u 6 (t)||L<-(n) < ||uo|U« ( n)e-^-* + - / M 2 (T)\\w £ (T)\\ L o C{n) e-— {t ~ T) d T . 

£ Jo 

Thus, using inequality (A.ll) and plugging the above estimate into (A.14) yields 



C w (t) <^ + ( HIU=o (n) - x J e~ d9t + <%o||L~(n) \ e-^-^e'^dr 
(A. 16) dg V t dgJ Jo 



d h -\-d 
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Changing the order of integration, we can simplify the last term on the right-hand side 
Jo 

1 4 w e(0\\L~ { n) c _ dJt _ e)cia _ dc) z dc _ f* £\\w e (Q\\ Loo{ n) c(a ^ c ±d£ ( ^ t)dc 
db + d — ed g J Q db + d — ed g 

hence, for < e < ^j-^, using inequalities (A. 12) we obtain 

(A.17) lMt)lUcc (n) <C + ^" L " (n j [\^Hwm\^ i u ) e- d ^dC, 

db + a — £d g Jo 

where C = C(||u>o||L°°(n), ||^o|U°°(fi)) is independent of T and of e. Finally, the Gronwall 
inequality applied to (A.17) implies the estimate 

(A.18) IK(0IU°°(fi) < Ce Cl{T) for all 0<t<T 

and for all sufficiently small s > (e.g. e G (0, (<4 + d) / (2d g ))) ) where positive constants 
C and C\(T) are independent of e. 

Finally, estimate (A.18) applied to inequality (A. 15) implies an analogous bound for v £ . 

□ 

Theorem A. 2. Let (u e ,v £ ,w £ ) be a solution of system (A.l)- (A. 5) with sufficiently small 
e > (e.g. e G (0, (<4 + d)/(2d g ))) and (u,w) be a solution of the corresponding system 
(A.6)-(A.10). For each T > 0, there exists a constant C(T) > independent of e such 
that 

(A.19) max \\u £ (t) - u(t)\\ LOO{n) < C(T)e, 

(A.20) max \\w £ (t) - iy(f)||i» (n) < C(T)e. 

Additionally, we also have 

(A.21) / \\v £ (t)-v(t)\\ Laa{n) <C(T)e, 

Jo 

where v is given by equation (A. 6). 

Proof. Letting a = u £ — u, (3 = v £ — v and 5 = w £ — w, we obtain by the Taylor expansion 
the following system 

(A.22) d t a = f(u £ , v £ ) - /(it, v) = -d c a + f x a + f 2 /3, 

(A.23) d t e(3 = g(u £ , v £ , w £ ) - g(u, v, w) - ed t v = -(d + d b )/3 + g x a + g 2 5 - ed t v, 

(A.24) d t 8 - DA5 XX + (d g + g 2 )S = h x a + d(3, 

supplemented with the initial conditions 

a(x,0)=0, 5(x,0) = 0, (3(x,0) = v (x) -v (x) 

with t>o obtained from uq and wq via formula (A. 6), and with the Neumann boundary 
condition for 8(x, t). In equations (A.22)-(A.24) the following coefficients 

f _ df (av) 2 df (au) 2 



du (u + v) 2 ' dv (u + v) 2 ' 

dg dg dh 2 dh 

g 1 = — = 2uw, g 2 = — = -— = u , h x = — = -2uw 

ou ow ow au 
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are calculated in certain intermediate points and are bounded independently of e due to 
Lemma A.l. 

The proof is divided into three steps. 

Step 1: First, applying the comparison principle to the parabolic equation (A. 3) with 
the Neumann boundary condition and with the zero initial datum we obtain the estimate 

||5(-,*)lk~(n) < C s (t) for every t G [0,T], 

where C$ is a solution of the Cauchy problem 

C s (0) = 0, 

— Cs + (d g + ||#2|U°°(nx[o,T]))C<5 = ||^i|U°°(Qx[o,T])||«(i)IU°°(n) + d||/3(£)||i/x>( fi ). 

Since 

C&{t) = /"%-(^+ll^ll^c«x [ o, r]) )(*-0(|| / , 1 ||^ (0x[o:r]) || ct(t) || ioo{n) + ^11/5(^11^^)^7- 
Jo 

using the Young inequality for a convolution and the estimate e _ ( ds+ " 92 " L ° ( 0x [°< T i)^* _ ' r * ) < 1, 
we obtain 

(A.25) ||5(-,i)IU°°(ft) ^ ||C5lU°°(o,t) < 1 1 'ii 1 1 (nx [o,Ti) 1 1 « 1 1 x, 1 (o,*;i«> (fi)) + d\\fi |Ui(o,t ; i°°(n)). 
Step 2: Next, we estimate the solution (3 = f3(x,t) of equation (A. 23). First, note that 

¥- {d+db) "\W m <-7^-r for all t > 0. 
d + db 

The solution of equation (A. 23) satisfies the formula 

P(x, t) = p(x, Q)e-^ + J e-^ ( *- r) \ -d T v + -( gi a + g 2 S)J dr. 

Consequently, the Young inequality yields 

||/3|U 1 (0,t;L oo (fi)) 

< (11/5(0) IU°°(fi) + ||^IUi(o,*;i°°(n))) C £ 
(A. 26) + C||5'i||L O o( nx [ 0iT ])||a|| L i( 0ii;LO o(Q)) + C 1 1 ^2 1 1 (rix [o,? 1 ]) 1 1 ^ 1 1 i-i (o,«;z,°° (n)) 

<Ce + C (||<7l||z,°°(nx[0,T]) + t\\hl\\L°°(nxlO,T])\\g2\\L°°(nx{0,T})) || Orjj ii(0,t;X,°° (fi)) 

+ 1 1 ^2 1 1 (n x [o,r] ) 1 1 ^ 1 1 i 1 (o,* ; i°° (n» 5 

where the last inequality results from (A.25). 

Step 3: Finally, we estimate the solution a = a(x,t) of equation (A. 22). Note that 
a(x, 0) = and 

a(t,x) = ~ J (/2(r)/3(r)e- dc( ^ r)+ ^ /l(?)d «) dr. 

Thus, using the Young inequality again we obtain the estimate 

(A.27) ||a|U°°((o,t)xn) < Ce( a ~ dc > T \\f 2 \\L°°(n)\\P\\L 1 (o,t;L°°(n)) 

as well as 

(A. 28) ||a|U 1 (°>*; iOG ( n )) — ^ll/^IU 1 (0.*;- E ' oc (O))- 

Inserting inequality (A. 28) into (A. 26) leads to 

(A.29) ||/3|Ui(o,tii<»(n)) < Ce + Ct\\f3\\ L i mLaam . 
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For t < t = 2Q we conclude that 

(A.30) ||/3|Ui(o,tiL»(n)) < Ce for all t < t . 

Since every t £ (to^] can be reached after a finite number of steps, estimate (A.30) 
holds for every t £ [0, T]. Furthermore, inequality (A.30) applied in (A. 25) and (A. 27) 
completes the proof of estimates (A.19)-(A.21). □ 

Remark A. 3. To obtain a better estimate of /3, one should construct an initial value layer, 
since (3\ t=0 ^ 0. 

Appendix B. Constant steady states 

Here, we consider the space homogeneous solutions of the two-equation model (A. 6)- 
(A.10) which satisfy the corresponding kinetic system 

(B.l) u t = ( " W d c )u, 

(B.2) w t = —d q w — , ^ b u 2 w + k , 

d b + a 

where a, d c , db, d, d g , k are positive constants and we have always assumed that a > d c . 
The structure of constant steady states of this system is the same as of the original three- 
equation model and can be characterized by the following lemma (for the proof see [14]). 

Lemma B.l. Let O = 4gL ( dc - | d b (d b + d). If nl > Q , then system (B.1)-(B.2) 

a \a - d c J 

has two positive steady states (u_,w_) and (u + ,w + ) with 

(B.3) u± = -{db + d) and w± = — . 

a — d c w± 2d g 

Theorem B.2. Let (u-,W-) and (u + ,w + ) be positive steady states of system (B.1)-(B.2) 
given by (B.3). Then (u + ,w + ) is always unstable. While (u-,W-) is stable, except for the 
case 

— {a-d c )-d g >0, -<1 and k < _ 

where 

— (a- d c ) 

P=dT >h 

— (a - d c ) - d g 

Proof. Let (u, w) denote a steady state of (B.1)-(B.2). From direct calculations, the 
Jacobian matrix J at (u, w) of the nonlinear mapping defined by the right-hand side of 
(B.1)-(B.2) is of the form 

/ dc, , x (a ~ d c f _ 2 \ 
a a{db + a) 

„ dbd c db _ 2 

v ~ 2 — t ~ d 9 ~ -rv~A u I 

\ a — d c db + a / 

We know (cf. Remark 5.2) that 



J 



UNSTABLE PATTERNS IN REACTION-DIFFUSION-ODE SYSTEMS 



27 



(i) if 



(B.4) <k( a -d c )-d g --^-u 2 <0 and - ^ (a - d c ) + 4 ^ f c) u 2 > 0, 



then all eigenvalues of J have negative real parts; 
fii) if 



(B.5) 



[a - d r ) - d n 



d 



b ,-.2 



d b + d 



u > or 



d a d, 



a — d c ) + 



d b d c (a - d c ) 



u 2 < 0, 



a(d b + d) 

then there is an eigenvalue of J which has a positive real part. 

Step 1. First, we show the stability of (u + ,w + ). Using estimates (B.3), the second 
inequality of (B.5) can be written in the form 

ra «n d b (d b + d) ( a - d c \ 2 2 

(B - 6) — 5— (— 1 <w 

Note that the left-hand side of (B.6) satisfies 



d b (d b + d) ( a — d c \ 2 



d n 



d r 



Adf 



and w + satisfies (K /(2d g )) 2 < iwi < (no/d g ) 2 . Therefore, it follows from the assumption 
K 2 Q > that 

Idf < idf < W+: 

which implies that the steady state (u + ,w + ) is unstable. 

Step 2. Next, we show stability of («_, The second inequality of (B.4) is equivalent 



to 0/(4dg) > w 2 . The latter inequality holds true, since w 2 _ 



and 



- f 2«5 - 2Koy«o - © - ©J = - - y'/t 2 , - 0J > 0. 

Using (B.3) and the relationship d b (d b + d) ^7rz^~j = jf-, the first condition of (B.4) 
becomes 



w_ 



— (a - d r ) - d n 



< 





Ad g 



If —(a — d c ) — d g < 0, then the above inequality always holds. 
Assume —(a — d c ) — d g > 0. Note that w 2 _ is given by 



w 



2k 2 -2/€o V / ^o- - 



4^ 



Therefore, it is sufficient to show the following inequality 

d, 



(B.7) 



2«g - 2hZ JrZ 2 Q -0-0 



— (a - d c ) - d s 
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The left-hand side of (B.7) is 



2Kt 



Hence, if 
(B.8) 



2k \/ Kg 



(a - d. 



— (a - d c 
a 



0-0 



— {a - d r ) - d a 



2Kn 



2k \. //d-0-0 



d„ 



2 Avn 



2kq\/ Kq 







2n r 



2k \ /k§ -0-0 



2k c 



2k \ kZ - 



d g Q. 



< 0, 



then inequality (B.7) holds true. Noting ^(a — d c ) — d g > 0, we obtain, from (B.8), that 
(B.9) 



what is equivalent to 
(B.10) 



p 



< 



— (a — a c 
a 



— (a — d r 



-0 = /30, 



d,, 



-0 < «oV«o 



0. 



If (3/2 > 1 and < Kq < (/3/2)0, then inequality (B.10) is always satisfied since the right- 
hand side of (B.10) is positive. The remaining cases are (i) /3/2 < 1 and (ii) j3/2 > 1 and 
> (/3/2)0. In the cases (i) and (ii), the both-sides of (B.10) are positive. Therefore, 
we calculate the square of both sides of (B.10) and obtain 

CB.li) i(h) e<Kl 

If $ > 2, then 0/2 > /5 2 /(4(/3 - 1)), while (3/2 < (3 2 /(4((3 - 1)) if /3 < 2. Therefore, the 
inequality (B.10) holds in case (ii). In case (i), (B.10) is satisfied under the condition 
(B.ll). □ 
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